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Abstract 



In principle every excitation acquires a finite lifetime in a hot system. This 
nonzero spectral width is calculated self-consistently for massive fermions cou- 
pled to massless scalar, vector and pseudoscalar bosons. It is shown that the 
self-consistent summation of the corresponding Fock diagram for fermions 
eliminates all infrared divergences although the bosons are not screened at 
all. Our solutions for the fermion damping rate are analytical in the coupling 
constant, but not analytical in the temperature parameter around T = 0. 
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I. INTRODUCTION 



Recently, the problem of the damping rate of a fermion moving through a hot system, 
i.e., the width of its spectral function, has attracted much attention This damping rate 

determines the radiation spectrum |]10| as well as the non-equilibrium transport properties 



of a relativistic plasma |[TT[], hence is an interesting quantity for a wide range of physical 
systems: Relativistic heavy-ion collisions, fusion reactors and the early universe are some 
examples. 

The damping rate of fermions interacting with massless bosons is an especially interesting 
problem in this framework. From the phenomenological point of view this importance is due 
to the fact that all fundamental forces of nature involve massless gauge bosons. From the 
theoretical standpoint a similarly big interest is due to the calculational difficulty associated 
with this particular problem. 

This difficulty is an infrared divergence, i.e., if interacting with strictly massless bosons 
a fermion may emit an infinite number of such bosons with infinitely small energy. Such 
a concept is of course meaningless in the real world, and consequently the solution of the 
problem has to connect the abstract world of quantum physics with reality. We may in- 
deed consider the solution as almost crucial for our understanding of nature by means of 
theoretical modeling. 

Such an infrared divergence occurs already at zero temperature in QED as well as in 
QCD. The connection to reality we mentioned, and which is conventionally used to remove 
this unphysical feature, is the finite energy resolution of every physical measuring apparatus. 

At nonzero temperature however, the question of the gauge- in dependent, causality pre- 
serving and therefore physical removal of the infrared divergence is far from settled. With the 
present paper we wish to contribute a new answer to this question, based on a perturbative 
expansion in terms of unstable asymptotic states. 

Our contribution is complementary to the most common method to remove infrared di- 



vergences at nonzero temperature, i.e., to the method of hard thermal loops |T2[. In this 
framework a partial resummation of diagrams contributes a small "thermal" mass to parti- 
cles, whereas their thermal scattering is usually neglected. However, this thermal scattering, 
or Brownian motion, is the fundamental difference between a system at nonzero tempera- 
ture and a system at T = 0, since it leads to a finite lifetime of any excitation. Indeed, one 
of the few mathematically rigorous theorems of finite temperature quantum theory, due to 



Narnhofer and Thirring, has established this once and for all |13|: "Particles" with infinite 



lifetime can exist in a hot system only if they do not interact. The use of particle- like excita- 
tions in a perturbative scheme, which are asymptotically stable and only have a temperature 
dependent "mass" is therefore unjustified a priori and can only be obtained as a limiting 
case at the very end of all calculations. 

The deeper mathematical reason is that the symmetry group of space and time in thermal 
states is not the Poincare group as in the vacuum state, but rather a product of SO (3) and the 
four- dimensional translation group . Its irreducible representations can be interpreted as 
having a continuous mass spectrum, i.e., the concept of a mass-shell (and thereby of stable 
asymptotic states) is not well-defined at nonzero temperature. 

The solution of this problem has been pointed out by Landsman [|I5| , using earlier work of 
Licht [T^ and Wightman [|l^ and leading to a perturbative expansion in terms of generalized 
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free fields without a mass-shell. To account for the proper temporal boundary conditions 
in a thermal system, these generalized free fields have to be embedded in a description with 
doubled Hilbert space, i.e., the single-particle propagators are 2x2 matrices [1^. Two 



flavors exist for such a formalism, the Schwinger-Keldysh (or closed time-path) method ||T9 
and the method called thermo field dynamics (TFD) |2 



Within the latter method this merger of different aspects of finite temperature field 
theory has been discussed in a recent overview [^, we will follow this description closely 
with our notation. Although in general we prefer TFD over the closed time-path method 
due to its technical simplicity, the present paper is formulated in a way that is independent 
of the choice of the formalism. 

Within this framework, we now consider the problem of the fermion damping rate. The 
paper is organized as follows: First, we make some general remarks on propagators in terms 
of their spectral function. These remarks are then used to obtain well-defined approximations 
to the spectral functions of bosons and fermions in a hot medium. 

In section III we briefly discuss the problem of calculating one-loop diagrams with ef- 
fective propagators, and how they may be used for a self-consistent determination of the 
spectral function. In the following we drop some of the previous results, concentrating on 
the fermion spectral width only: Sections |V|and|V| are devoted to analytical approximations 
for the Fock self energy function. 

In section ^ we assemble these results, and consequently obtain approximate solutions 
for the self-consistent Fock approximation, i.e., for the full fermion propagator calculated 
with a one-loop self energy function at finite temperature. Finally these results, which are 
different from those obtained by other methods, are used to draw some conclusions. 



II. SPECTRAL FUNCTIONS 

As appropriate for a description of dynamical phenomena in a thermal system, we em- 
ploy a formalism with 2x2 matrix valued Green functions. In thermo field dynamics (TFD) 
this "doubling" arises naturally because TFD contains two different (anti-) commuting rep- 
resentations of the canonical (anti-) commutation relations. 

In the Schwinger-Keldysh method the doubling is introduced because a full description of 
statistical systems can be achieved only when using causal and anti-causal Green functions 
as well as Wigner functions. Consequently, the full propagator matrix is 



where (■) denotes the statistical average, T[-] the time ordered and T[-] the "anti-time or- 



dered" product. The same matrix valued propagator is obtained in TFD (see pl|j2^ for 
detailed discussions of this equivalence). 

By construction the propagator matrix obeys the linear relation S^^ + S^^ — 5^^ — = 0. 
Furthermore, in an equilibrium state it has to fulfill the Kubo-Martin-Schwinger condition 



23| , i.e., an anti-periodic boundary condition in the imaginary time direction. We will 



henceforth formulate all our considerations in momentum space, the KMS condition then 



reads [18,21,22 
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{l-nF{po))S'\po,p)+nF{po)S^\po,p) = 



(2) 



S^"^ and S*^^ are the Green functions without time ordering (Wigner functions) and np{E) 
is the fermion equihbrium distribution function at a given temperature, the Fermi-Dirac 
function 

ME) = . (3) 

An almost identical KMS relation holds for bosons, but here the fields are periodic in the 
imaginary time direction. Consequently, occupation number factors are given as Bose- 
Einstein function 

nsiE) = -^J—^ . (4) 

Henceforth for simplicity we consider only the case of zero chemical potential. 

In interacting systems, the KMS boundary conditions lead to a causal propagator which 
has a cut along the real energy axis. Its analytical structure therefore is not easily under- 
stood, and especially when combining several of such propagators in a perturbative scheme 
one has to implement more or less complicated cutting rules to understand pieces of diagrams 
physically @,|2|]. 

Although this may be a matter of personal taste, we therefore deem it the safest method 
to use only retarded and advanced propagators, and to condense the matrix structure of eq. 
(HI) into the vertices where such propagators join. This amounts to a diagonalization scheme 
for the matrix (p, which has been described in several publications pT| , P^ . 



A. Retarded and advanced propagators 



The retarded and advanced propagators are, by definition, analytical functions of the 
energy parameter in the upper or lower complex energy half plane. Analytical functions 
obey the Kramers-Kronig relation, and this implies that the retarded propagator is known 
completely if only its imaginary part (or spectral function) is known along the real axis. 
Hence, for arbitrary complex E 

oo 

S''"'{E,p)= jdE' Af{E\p)^-^^^. (5) 

— oo 

Trivially, the spectral function is recovered for real 

Af{E,p) = T-Im(5^'^(i?,p)) = ^ {S^iE,p)-S\E,p)) . (6) 

We consider the diagrammatic rules to combine these propagators in the calculation of phys- 
ical quantities to be well established. For completeness we summarize the diagonalization 
of eq. (|1]) in terms of the spectral function Af{E,p) as 
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S^'''\po,p)= JdEAp{E,p) X 

— oo 

/ 1 \ 



Po- E + ie 



B{np{E)) . (7) 



V po-E-ie / 

Ts is the diagonal Pauli matrix, and the transformation matrix B is 

B(nAE)) = ( (/ - "-'^)' -"-'^) ) , (8) 

It is a matter of a few hnes to show that this matrix-valued propagator obeys relation 
for each spectral function. 

A similar relation holds for the boson case, see |PT| , |22| for details of the corresponding B. 
For both cases, the transformation matrices B have a special meaning in the TFD formalism, 
where they play the role of a thermal Bogoliubov transformation. However, the explicit form 
(^) of the propagator is the same in the Schwinger-Keldysh method. 

B. Normalization and Locality 

The spectral function has two features which are intimately related to fundamental re- 
quirements of quantum field theory. We will comment on these only for the fermionic case, 
and refer to for the case of bosons. Firstly, the quantization rules for fields (which we 



will use in the free as well as in the interacting case), 

{^{t,x),ij\t,y)}=6'{x-y) (9) 
require that the spectral function is normalized 

oo 

JdETT[-f''AF{E,p)] = 2. (10) 



The second important feature of the spectral function is that its four dimensional Fourier 
transform into coordinate space must vanish for space-like arguments. This is equivalent 
to the Wightman axiom of locality, i.e., field operators must (anti-)commute for space-like 
separations in Minkowski space 



if X — y space-like . (11) 
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In an interacting many-body system, we may very well expect non-locality in a causal sense: 
Wiggling the system at one side will certainly influence the other side after some time. 
The locality axiom ensures that this influence does not occur over space-like separations, 
i.e., faster than a physical signal can propagate. Thus, to distinguish between the causal 
non-locality and the violation of the locality axiom, we will henceforth denote the latter a 
violation of causality. 

It is trivial to establish that normalization and locality axiom are satisfied for free 
fermions of mass M. Defining their on-shell energy as ujp = a/p^ + , the free spectral 
function is obtained as 

A'piE, p) = + p7 + M) sign(E) 6{E^ - ^J) . (12) 

Without loss of generality we set y = in the free fermionic anti-commutator function, and 
obtain 

C-o (x, 0) = {id^-f^ + M) Z{xo, x) where 

Z{xo,x) = ^ (^6{xl - x') - Q{xl - ^")^^y=^, MM^xl-x^)^ , (13) 

where Ji is a Bessel function of the first kind. Clearly this is zero for space-like arguments, 
i.e., for \x\ > \xo\. However, from this prescription follows that for a general spectral 
function the locality axiom is not automatically guaranteed, a careful check is necessary in 
any application involving spectral functions. 



C. Ghost poles 



For an interacting system, the full fermion propagator and thus also the spectral function 
is defined in terms of a self energy function S: 



=1=— Im 

71 



iP0,P)) ^ 



(14) 



Note, that also the self energy function is Dirac matrix valued, i.e., this equation contains 
the inversion of a 4 x 4 matrix. 

However, the real part of the full propagator is not given by the real part of the quantity 
in square brackets of (14). The reason is, that according to Weinberg's theorem the self 
energy function in a relativistic theory behaves as 



S(po,p) ocpo (log (1^)) if 



Po 



00 



(15) 



i.e., it is more than linearly divergent. This implies, in general, the presence of unwanted 
poles. For fermions, these Landau ghost poles appear at four points in the complex energy 



plane |2^, p. 636]. 

Since their appearance contradicts the definition of the retarded/advanced propagator 
as being free of poles at least in one half plane, they must be considered unphysical. For 
sufficiently large coupling, they have a big influence on the propagator along the real axis. 
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while for small coupling they are far away from the physical region. Hence, their consistent 
removal is mandatory for strongly interacting systems, and the only way to do so is by 
calculating the real part of the propagator by a dispersion integral according to (^). 



It was pointed out in ref. |30|, that this indeed conforms to a non-perturbative correction 
to the propagator, in the sense that the correction term is not an analytical function of the 
coupling constant. In general, it will be of order exp{—l/g'^) instead. Hence, eqn. (|^) is 
the only safe method to calculate a relativistic propagator, and the input needed for this 
calculation is a physically meaningful spectral function which obeys all the rules discussed 
here. 

To demonstrate briefly how the removal of unphysical poles is achieved through a dis- 
persion integral, we discuss a bosonic toy model spectral function of Lorentzian type which 
has a peak along the real axis at some value Up 

Ao{E,p) = l \ (16) 

By definition, this spectral function may be used only for real energy arguments, and does 
not take on negative values. If this were violated, the probabilistic interpretation of the 
quantum (field) theory behind this formulation would be invalid. 

The above function has a trivial expansion into simple poles in the complex plane, 

ME.P) = ^ U ^— - ^ ) . (17) 

We therefore calculate the real part of the corresponding propagator through 

G^{po,p)= [dEAo{E,p) (18) 

—00 

The integrand is expanded into simple poles, 
1 1 



E - ujp±i'jp po - E + ie 

I ( I + I ^ . (19) 

Po - (ujp T ijp) \E - Up ± I'jp po-E + ieJ 

The i?- integration is convergent, and closing the integral in the upper complex energy half 
plane then picks up the residues of the corresponding poles. First consider the case Im(po) > 
0, where the equal sign is permitted because of the ie terms in the denominator of the 
dispersion integral. Clearly, the integrand then has three poles in the upper complex half 
plane, but two of the residues cancel. 

In the second case Im(po) < the integrand has only one pole in the upper complex half 
plane. Assembled together the integral gives 

^ if Im(j9o) > 



G^iPo,p) = <'''~^\'~'^'^ , (20) 



Po - {Ep + i7p) 



if Im(po) < 
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which is free of poles on both sides of the real axis. A similar effect occurs, when a self energy 
function according to (^) is used, i.e., when ghost poles would be present in a perturbative 
propagator: The dispersion integral always shifts complex (ghost) poles to the unphysical 
Riemann sheet. 

As a matter of fact it has been shown very long ago in quantum field theory, that this is 
the only scenario compatible with causality: The retarded as well as the advanced propagator 
must be free of poles in the complex energy plane on both sides of the real energy axis pi 



D. Approximate spectral functions 

The toy model spectral function of eqn. ([T6|) cannot be used in a relativistic calculation, 
since it propagates only states of positive energy. More useful is a symmetric ansatz, which 
for simplicity we first discuss for the bosonic case: 

Ab{E, fc) = (A(E, k) - A(-^, k)) 

= 1 (21) 

with Vl\ = ul + 'jI. This spectral function is normalized according to the canonical field 
commutation relations, i.e., 

oo 

jdEEAB{E,k) = ^. (22) 



Note, that in this ansatz the dependence of 7^ on the momentum of the relativistic boson 
is completely arbitrary, and one may also introduce a more general relationship between uuk 
and momentum than in the free case. 

It is a matter of a few pages to obtain the full propagator from the dispersion integral 
over this spectral function as 

^""(^o, k) = ^ (- , , - „ , \ , , ] \og{-ko - ie) 



+ 



2niUk V (^0 - ^^kY + ll ih + uJkY + ll 
1 / log ( -tUfc + i7fc ) log ( - i7fc 



AiTiUk \ ko-uJk + i'jk ko-LOk- i7fc 
1 / \og{uJk + i7fc) log(u;fe - i7fc 



ATTiujk \kQ + ijJk + i7fc k^ + Uk- i7fe 



(23) 



This propagator is valid everywhere in the complex plane, and it is by construction free 
of poles on the physical Riemann sheet. On the real axis, it has the imaginary part 
—TiABikoi k), hence reproduces the ansatz for the spectral function. 

The question is now, how one may relate such an approximate spectral function to the 
general full propagator. Clearly it amounts to the approximation of the denominator of a 
relativistic retarded boson propagator along the real axis as 

k^k'' -m^- n^(fc) = {ko - {iOk - hk)) {ko + {uJk + hk)) , (24) 
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where 11 is the boson self energy function, and solving this equation for ujk, 7fc. An example 



for such an approximation has been given in ref . |2T[ , and there the momentum dependence 
of as well as of 7^ was quite strong for the case of a pseudoscalar coupling between bosons 
and fermions. 

Let us remark at this point, that a general momentum dependence of the 7fc-parameter 
does not guarantee the causality of the model, i.e., it may violate the locality axiom as 
discussed above. For a physically correct model, it is mandatory that the Fourier transform 
of the spectral function vanishes for space-like coordinate arguments. Indeed one may show, 
that for the ansatz eq. this is achieved when considering an ujI which is quadratic in k 
and a constant 7^. 

This ansatz may be extended to the fermionic case. The spectral function then is slightly 
more complicated because of its Dirac matrix structure. In eqn. (|10|) we already specified 
the normalization for a fermion spectral function. The simplest ansatz compatible with this 
normalization, as well as with known properties of fermion systems, would be to multiply 
the spectral function from eq. (pT]) by a factor (-£^7° + ujp). 

However, as one may show, such an ansatz (used in an earlier version of the present 
paper) violates the locality axiom, therefore is not useful for our purpose. We therefore 
resort to a more complicated spectral function, 

7, 7° {E' + ujI + 7,^) + + 2i^M 

A.F\i^iP) = 

ujpl^ +PJ + M -cjp7° + PJ + M 



4:7iiUp \ E — Up — i'-yp E + Up — i'jp 

Wp7° +P7 + M ^ -Up-f^ + p-f + ^25) 



E - LOp + i'jp E + ujp + i7p 

where in principle 7^ may be momentum dependent, fip^ = cUp"^ + 7^^ and M is a constant 
which may be different from the physical fermion mass in the vacuum state. 

We now again turn to the question of consistency. As argued in the previous subsection, 
the Fourier transform of the spectral function has an important physical interpretation. We 
observe that in the special case where 7 is independent of momentum the 4-dimensional 
Fourier transform of eq. (pS]) is 



Cpix, 0) = e-^ l^ol (^h^^ + Z{xo, x) , (26) 

with the function Z as defined in eq. (0). Hence, for momentum independent 7 and M, 
the locality axiom is fulfilled and interactions in the system are causal. We therefore limit 
the use of the spectral function (|25| ) to this approximation, i.e., 

7p =7 = const. 

=^ uj{p + kf= + 2\p\\k\T] + , (27) 
where rj is the cosine of the angle between the two momenta p and k. 
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III. FERMION DAMPING RATE 



The next step is to use spectral functions, together with the definition of the complete 
propagator by dispersion integral, in a "perturbative" expansion at finite temperature. We 
find, that "perturbative" is certainly not the correct label for such a skeleton expansion, 
because in general due to the properties of the dispersion integral the results will not be 
expandable into power series in the coupling parameters (see section II. C). Indeed up to 
correlation diagrams and vertex corrections, a perturbative expansion in terms of full prop- 
agators is exact already at the one- loop level [^, pp. 476]. 

For simplicity, we begin the discussion with a simple scalar coupling of a boson field 
having the spectral function Ab to a fermion field having the spectral function Ap- In eqns. 
(PPP and (^Op, some generalizations are investigated. 

The "one-loop" diagram for the fermion self energy, the Fock diagram, is given by the 
integral 

S^(j>o,p) = 

CO 

— oo 

This function is split into real and imaginary part as = ReE — ivrF, and for the latter we 
obtain 

r(po,p)= (29) 

oo 

^7(2^ yc?^^F(^ + Po,P + fc)^B(^,fc) {nniE) + npiE + po)) . 

— CO 

To fulfill the Kubo-Martin-Schwinger boundary condition (Q) for the fermion propagator 
including the above self energy function, it is absolutely essential that both propagators in 
the above expression have the same equilibrium temperature. In other words, one may not 
disregard the fermion distribution function up^E + po) in the above expression in any case, 
as was done e.g. in ref. [Q]. This fact will play an important role in the next section as well 
as in sect. VLB. 

The corresponding "one-loop" result for the retarded scalar boson self energy is, accord- 
ing to ref. [pll] 



n^(A:o,fc)=/y"^ jdEdE' 

— oo 

^WE,,,U)ME',,)^{^^^lf^). (30, 

The imaginary part of 11^ = Rell — ma is explicitly 

oo 

-(^o,fc)=//^ jdE 

— CO 

Tr [AF{E,p + k)AF{E - ko,p)] ME - k^) - upiE)) . (31) 
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The above imaginary part vanishes at ko = 0, more generally when the external energy 
parameter equals the boson chemical potential. This relation is violated if different propa- 
gators for the fermions are used on the two legs of the diagram, e.g., when inserting different 
"thermal masses" for the two legs This violation in turn leads to a violation of the 



Kubo-Martin-Schwinger (KMS) boundary conditions for the full boson propagator, i.e., it 
is not an equilibrium Green function. 

In principle, one could now introduce a self-consistent calculation scheme: The expres- 
sions ( p8D and (|30D are calculated with an approximate spectral function, then used as input 
to eq. (|I^) and a similar equation for the full boson propagator in terms of its self energy. 
This gives new spectral functions, which may be used again to determine self energies. Such 
a scheme has been used in refs. [^,|^, and based on these papers one might hope that a 
few iterations are enough to determine the spectral function numerically quite well. 

However, as we have pointed out, there is no way to guarantee that in such a scheme 
causality (in the representation of the locality axiom) is preserved. Moreover, in an entirely 
numerical scheme it would be impossible to point out the path to the solution of the infrared 
problem. Two strategies exist which may be used to circumvent this difficulty. 

The first strategy is to enforce locality after each iterative step. This can be done by 
folding the resultant spectral functions with the Fourier transform of G(t^ — x^), which 
makes this scheme numerically much more difficult than even the iterative procedure used 



in 11,14 



The second strategy is to use a simple causality-preserving parametrization of the spectral 
functions. The parameters are then determined self-consistently in a closed scheme. Apart 
from the possibility to achieve analytical approximations, this strategy also has the virtue 
that one can use it as an input to the first one. Thus, for obvious reasons we will follow this 
second course. 



A. Approximations to the boson spectral function 

As was pointed out above, and is well known for many years, the boson propagator 
may have isolated poles only on the real energy axis, but not away from it on the physical 
Riemann sheet Furthermore it may exhibit cuts along the real energy axis due to self 
energy corrections with continuous nonzero imaginary part like in eq. (|3T|). 

In absence of condensation phenomena these self energy corrections have zero imaginary 
part when the boson energy parameter is equal to the chemical potential. Consequently, 
the product of a Bose-Einstein distribution function (^) and a continuous part of the boson 
spectral function is always infrared finite. 

The only possible source for an infrared divergence therefore are poles on the real axis, 
corresponding to particles with infinite lifetime. We have argued before, that in principle 
such poles does not exist in a finite temperature system |T^, and thus it is one on hand 



clear that infrared divergences in finite temperature field theory are not present in the full 
theory but a mere artifact of perturbation theory. 

On the other hand it is nevertheless instructive to study how these divergences are 
removed within a proper calculational scheme. We therefore chose the worst case for a 
boson spectral function, corresponding to a kind of electromagnetic interaction without 
screening or damping: 
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A'^ = -iSiE-k)-SiE + k)) (32) 

with k = \k\. It is obvious, that the product of this spectral function with the corresponding 
(particle plus antiparticle) Bose-Einstein distribution function is highly singular for — > 0. 

However, inserting this spectral function into the calculation of F gives at the lower 
boundary of the /c-integration 

/dk 
^ (47r 2T Af{po, p) + Oik')) , (33) 



i.e., the integral is infrared finite as long as the fermion spectral function is finite. That it 
is also ultraviolet finite is not obvious from the expansion made here, but follows from eq. 

We therefore conclude, that the use of a proper finite temperature fermion propagator, 
which obeys the Narnhofer-Thirring theorem and therefore has no isolated poles anywhere 
in the complex energy plane, is completely sufficient to remove all infrared divergences from 
the calculations. In particular, no screening or damping of the bosonic interactions is needed 
to achieve this effect. 

Introducing a modified (i.e., less singular) spectral function for the bosons of course 
modifies the result for F quantitatively, but only one direction for such a modification is 
possible: The calculation of F is dominated by the infrared sector. Any less singular spectral 
function than ( p2D therefore diminishes the value for F, because it leads to a "smearing" of 
the Bose pole with some distribution. Consequently, we may consider the values obtained 
in the following an upper bound on the full calculation. 



B. Approximations to the fermion spectral function 

We now set up a self-consistency scheme for the determination of the fermion spectral 
function. Suppose, we really would approximate the full spectral function according to eq. 
(|l^) by the simple parametrization of eq. (p5|). This then would amount to equate the 
denominators on the real energy axis as 

{p, - V,{p) + ivrF;(j9))' - (M + S{p) - i7rr{p)f 
= {Po - i^p - hp)) {Po + i^p + hp)) ■ (34) 

The functions of momentum appearing here are the components of the retarded self energy 
function, split into real and imaginary part according to 

S^(p) = Re(S(p)) -i7rF(p) 
Re{E{p)) = Sip) + V,ip)j^ 

r{p) = r{p) + r;{p)Y. (35) 

In the above expressions, we have also split the self energy functions in a Lorentz scalar and 
a Lorentz vector piece, henceforth simply abbreviated as scalar and vector part. 
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We have stated above, that in general this may give rise to a momentum dependent 7^ 
as well as to a non-quadratic dependence of ujp on p. Consequently, to check for the locality 
axiom, i.e., the preservation of causality in the interactions of the system, requires great 
numerical effort and one would lose all the virtues that come with a simple parametrization 
oi Af- Let us therefore study first, how the choice of spectral function ( pSf ) affects the 
imaginary part of the self energy (^91) . Clearly it has only two independent components, 
since its Lorentz scalar and vector part are 

r^(p) = r^(po,P,7) 
rr(p) = ^r^'(Po,P,7) 

r(p) = r^^(j9o,p,7). (36) 

In the following, we want to study a slow massive fermion on its effective mass-shell, i.e., 
for pq = ujp = ^Jp^ + M^. We therefore make the approximation that the real part of the 
self energy function may be absorbed into the mass parameter M. 

For the scalar coupling we consider here, the self-consistency equation obtained from eq. 
fBl) then reads 



75 



7rr^(^, p, 75) + 7rr^^(^, p, 75) (1 - ^) ^ ■ (37) 



According to our assumption, the left side does not depend on the momentum p. While we 
certainly cannot hope that the momentum dependence on the right side cancels completely, 
we may nevertheless expect that it possesses an expansion of the type 

75 = ^^o + ^^i^ + '^2^ + ... . (38) 

Consequently, our approximations are reasonable if the coefficients in this expansion are at 
least of the same magnitude - and they must be considered a failure when e.g. <^ |z/i|. 
To check this condition therefore will provide an a-posteriori test of our assumptions. 

It is quite simple to generalize this scheme to the exchange of a massless vector boson, 
if the boson propagator is taken in Feynman gauge p9|, pp.329]: The inclusion of Dirac 
matrices 7^, 7^ at the vertices of the Fock diagram simply multiplies by a factor -2 and 
F* by a factor 4. The resulting vector boson self consistency relation is 

7y = -2nT\u,p,^v)+^7iT'\uj,p,^v) f 1 + — • (39) 

It is worthwhile to note, that the ansatz of a momentum independent imaginary part of 
the fermion self energy ensures the gauge invariance of this self-consistent damping rate in 
the limit p — >■ 0. Gauge invariance even of the first momentum dependent correction term 
requires to introduce a vertex correction in eq. (^) to satisfy the correct Ward identity |^ 



Physically interesting is furthermore the result of massless pseudoscalar bosons, where 
one has 



7p = 7rT\uj, p, 7y) - ttT^\uj, p, 7p) (^1 + 



M , , 

- ■ («) 
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C. Angular integrations 



The integration over the angle between the two momenta k and p may be done analyti- 
cally. Using the above spectral functions, we obtain as the 7°-proportional imaginary part 
of the self energy 



oo 



47]-3 



X) 

Uk |fc7 (7° {{k + ujf -kuj + 7V2) h{k + u) 

+-i%k h{k + uo)) {nB{k) + npik + po)) 
- kj (7° {{k - + ku + 7V2) hik 



UJ 



+jyk hik - u) (uBi-k) + npi-k + po))^ . (41) 

/i and I2 are integrals over the angle between the momenta p and k, and we have abbreviated 
k = \k\, p = \p\ (see appendix A for technical details). 

A similar expression can be found for the second piece of the self energy function, 



r"{uJ,p,l) 

47]-3 



2 °° 

M-f- dk{k-f{k + uj) h{k + uj) {nB{k) + npik + po)) 







+ k-f {k - Lu) Ii{k - Lu) (uBi-k) + npi-k + Po)) } . (42) 



The infrared problem addressed above finds its representation in the fact, that in the limit 
7 — >^ 0, Ji and I2 behave like a negative power < —2 of k. This becomes obvious when 
expanding the integrand around the value k = 0. 

To check the conjecture of infrared finiteness made after eq. ([33|), we perform this 
expansion in the above expressions and obtain for the quantity in curly brackets of eqn. 



(PI), at the point po = u that we are specially interested in: 



Consequently, the infrared divergence is removed by starting with a small nonzero 7. 



IV. VACUUM STATE 

The crucial difference between the present paper and hard thermal loop calculations 
||l|-0 is, that we do not find it necessary to introduce momentum cutoff scales: The quantity 
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we are calculating is infrared and ultraviolet finite, and therefore no justification exists to 
neglect the "vacuum" contributions. These arise, because the Bose-Einstein and Fermi- 
Dirac distribution functions are nontrivial also for zero temperature, n^^{E) = —Q{—E) and 
= Q{-E). 

One may argue, that a "particle" in the vacuum generally is a well-defined state, i.e.. 



the NRT theorem |13| does not apply and therefore the damping rate should be zero. This 
is of course correct, but we are not using the ordinary vacuum propagators for the fermion 
field. Instead, we must calculate the contributions with a nonzero fermion damping rate as 
input, and only later may be able to find a self-consistent solution of zero damping rate. 

Using the above step functions in the expression for the imaginary part of the self energy 
leads to the simpler form 

r^vac(a;,p,7) = (44) 



oo 



jdkk-f {{{k - Luf + kuj + ^^2) h{k-u)+pkh{k-uj)) 60 



uj — k) . 
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This part of the self energy function has an asymptotic behavior according to Weinberg's 



theorem ||28|: The above integral is ultraviolet finite, but the corresponding real part is 
divergent and has to be renormalized. It is connected to the above imaginary part only 
through a subtracted dispersion relation. This asymptotic behavior of the self energy function 
for large values of implies, that the ghost problem discussed in sect. II. C becomes relevant. 

In the light of our introductory remarks to this section it is instructive to study the 
integral in case of zero 7 and momentum, where it becomes 

r..>o,o,o) = — -po pi^pi + M-) 0(bo|-M) 

- Po if Po 00 . (45) 
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This imaginary part vanishes at po = more generally at po = where we want to 
calculate the spectral width. We may put this into the simple form: Zero input 7 gives zero 
output V^^ for "on-shell" particles. 

However, this does no longer hold once we calculate the integral in (^) with a nonzero 
positive 7 as input: Then, we obtain a nonzero imaginary part also at the point Pq = uj. In 
figure 1, we have plotted r^^^(po, 0, 7) for two different values of 7: At pq = M, the imaginary 
part of the Fock self energy function is zero only in case 7 = 0, otherwise it is positive. 

One may therefore complete the statement above: Nonzero input 7 gives nonzero output 
r,,, even at po = ^• 

We now exploit the virtue of a simply parametrized spectral function for fermions, by 
expanding the continuous function F^ac around its value at 7 = 0. As we have seen above, 
this value is zero. However, it turns out that the imaginary part of the self energy function 
with pq = UJ is not an analytical function of 7: It has a cut in the complex 7-plane, starting 
at 7 = 0. Thus, a Taylor expansion in 7 is not possible, one has to perform an asymptotic 
expansion, explicitly taking into account the leading nonanalyticity. 

Various techniques exist for such an expansion, for the purpose of the present paper we 
isolated the leading terms by first substituting k ^ x'-f and then expanding the integrand in 
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powers of 7. It turns out, that for finite momentum p the scale for the expansion is u"^ — 
- which implies, that this asymptotic expansion is reliable only up to momenta p ^ M. 

Within this limitation, our result for the vector imaginary part of the Fock self energy 
function is 

r_(^,p,7) = -7-^<^ T-^— :o log H^— ^ -i + -iog' 



4tt^ I A uj"^ — p"^ \ \uj'^ — p^ J uj \uj — p 

7^ TT uj{uj'^ + 3p2) 



(tu^ — p^Y 



73 



24 (c;2 - p2)3 



1' 



c;2 — p2 



(2a;* + 18a; V + 12p*) log ^ } + 0(7*) , (46) 

where we have used p = \p\ and k = \k\. In figure |^, this expansion in 7 is compared to 
the full numerical calculation of eq. (P5|): With each additional order that is included, the 
quality of the approximation grows. 

An important aspect of this expansion is the occurrence of the quadratic term: In a 
naive view of the integral involved, it is not present because Ii{k — po) as well as l2{k — po) 
are odd functions of 7. However, as pointed out above, the function we are expanding is 
not analytical in 7 - and this effect causes the appearance of the quadratic contribution. To 
demonstrate this, figure 2 also contains a curve where this quadratic term has been left out. 
Obviously in this case the third order contribution makes the approximation worse instead 
of improving it. 

In view of the above results we conclude, that for not too large values of 7 the asymptotic 
expansion to second order is well under control. Henceforth, to keep the results simple, we 
therefore restrict the discussion to this next-to-leading log order. 

We now turn to the second piece of the imaginary part of the fermion self energy. Ac- 
cording to (|4^ ) one obtains in the vacuum 



r:i(a;, p, 7) = -m|-^ ldkk^{k~u)h{k~u)e{LO-k). (47) 



Performing the asymptotic expansion as above then gives 

rvac^.P>7 = 2 log ^ ^ +-log 

An-^ { 4: uj^ — p^ \ \oj^ — p'^ J p \oj — p 
7^71 M(3c;2+p2) 



(o;^ — p'^y 



7=^ M 

7 



a;(a;^ — p^Y 
-W{uj^ +p^) log 



Q (21a;* + 14c;V-3/) 



a;2 — p2 
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p 



We notice again, that as a sign of the non-analyticity in the parameter 7, a quadratic as 
well as a 7 log(7) contribution appear. 

V. THERMAL STATE 

We finally proceed by considering a thermal state, i.e., we calculate the full imaginary 
part of the self energy function according to (|29|) , which we split into the "vacuum part" 
calculated above and a temperature dependent part according to 

r(po, p, 7) = ^.Upo, p, 7) + ^Apop, 7) 

= r_(po, P, 7) + 4(po, P, 7)7° + r^^bo, P, 7) (1 + ^) ' (49) 

similar to the decomposition for the vacuum part only in eq. (pBI). For the moment it seems 
quite hopeless to obtain an analytical approximation to the full expression ( ^T] ) valid for all 
temperatures. However, we know that for not too high temperatures the Bose-Einstein and 
Fermi-Dirac distribution function may be expanded around their zero-temperature value. 



nB{k) + npipo + k) ^ exp{-k/T) [t + n+ ^piPo) + 77; ttt + ^piPof 



T 1 . . k f 1 
nsi-k) + npipo - k) ^ -exp(-fc/T) {^ + \ + 

+ 0(Po - k) (^-1 + np{po) + ^nF{po){l - npipo))^ • (50) 

The term with the —1 in the last line constitutes the "vacuum" part as obtained in the 
previous section. The other terms are ordered according to their dominance, i.e., the contri- 
bution proportional to T/k is the largest due to the strongly peaked nature of the function 
exp(— fc/T). The terms including the fermionic distribution functions are negligible for 
Po>T. 

It was stated in eqns. (|33| ) and (|i3|), that the introduction of a small nonzero 7 removes 
the infrared divergence from the integral. Thus, similar to the vacuum case, one obtains a 
result for Ft which cannot be expanded in a Taylor series around the point 7 = 0. 

We have described above how to obtain an asymptotic expansion around this value. Here 
this is again achieved by rescaling A; — > X7 in the integrand, and then expanding everything 
but the exponential factor exp(— X7/T) in powers of 7. 

To improve the calculation, yet another trick is introduced: First, the momentum integral 
is calculated using only the leading term T/k = l/{j3k) of the above expansion. The result 
is then improved by acting on it with a differential operator 

ri'f ^ ri( ^ ^ I rj.(a;,p,7) 

r,(o;,p,7) = r,(.;,p,7)-^-^-^^J T ■ 
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For both pieces of the self energy function this is a quite laborious task, for our calculation 
we heavily relied on symbolic computation using Mathematica*™. For the first piece one 
obtains 



rT(-,P,7) = ^|^iog^^J 



7 ^2 / / 2^2 X ^ (uj + p\ 17 20T2 



4^72—^2 \ \T'^{(jj'^ ~ p^) J p \uj — p J 3 3a;^ 



24T (Cj2-p2)2 



(^2_p2)2\^ 8j9 V V^(^ + p)y 12 



8p V V^(^ -P)) 12 

+ (3-^ + /)(^-^r^))}+Oh*), (52) 

where Cr = 0.57721... is Euler's constant. A comparison of this approximation to different 
orders in 7 with the full numerical calculation of the temperature dependent contribution 
is presented in fig. |[ As before, it turns out that the results are already quite reliable in 
second order of 7, provided the temperature is not too high. 

In case of the vacuum parts of the self energy function we were able to gather all terms 
of a given order in 7 in this asymptotic expansion. For the present temperature dependent 
pieces this is not possible, because of the above expansion: It automatically counts T and 7 
to be of the same order, hence is a simultaneous expansion in these parameters. As a small 
reminder of this fact appears a term 7T2 in the above expansion, and the 0(7^) must be 
replaced by a (9(7"^, 7T^) 

We find, that in second order 7 the asymptotic expansion is reliable up to 7 ~ T/2, 
whereas it is good up to the temperature in third order. In anticipation of the following 
section, we emphasize that the quality of the approximation grows tremendously for very 
small temperatures T ^ O.IM 

With similar quality we may also obtain the second piece, asymptotically given by 

A ^ flog i^,;^^ + ^ log f ^1 + 2C, 



A u'^ — p'^ \ \T'^{uj'^ — p^) J p \uj — p J 3 
-i'^TiMuj + 12T2 



24T (Cj2 _ p2)2 



log X +Cr- 



cu2-p2)3\^ 8p V \T{^^V)) 12 



8p V V^(^-P) / 12, 
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To this expression one may apply the same considerations as before: Although there is no 
piece oc 7T^, we may trust this expansion in second order 7 only up to 7 ~ T/2. 

The claims we lay on the accuracy of our asymptotic expansions are supported by fig. 
^, where we compare them with a completely independent numerical calculation We have 
checked the accuracy of our results with various values for momentum and temperature, and 
found it to persist up to momenta p ~ M. 

To be completely sure of our findings, we furthermore will drop the corrections of orders 
7'^: Unusual as it may seem for quantum field theory, we therefore have an approximate 
method with controlled accuracy in the next-to-leading-log order. 

VI. FERMION DAMPING RATE II 

Finally, we assemble the expressions for the self energy function that we have obtained 
in the previous two sections. Combining ( PB| ) with ( |5^ and (^Sp with (|55D, we then observe 
that all the terms of order 7 log(7) cancel ! 

This occurs for scalar and vector part of the self energy function independently, hence 
hold for all types of massless bosons that were considered in the equations (|37D - (^0]). We 
will comment on this in more detail below. 

As a consequence of this fact, the self-consistent damping rate for the fermion moving 
slowly through a hot medium, as defined in the equations (|37|) - (^Of), is the solution of an 
algebraic equation of the form 

Mco,P,T)-(^Mu;,p,T) 7 + /2Kp,T)72 = 0. (54) 

This equation has always two solutions for 7, but in general one of them is negative or very 
large. Hence it is obvious how to choose the physical solution of the above equation, and in 
the following we are referring to this one only. 

Our results do not exclude the possibility of corrections of order 7^ log(7) to this equation. 
The coupling constant appears only at one point, the three functions /o - /2 do not depend 
on it. /o is proportional to the temperature T, fi contains terms which are not analytical 
around zero temperature, i.e., terms of order T log(T). 

The solution 7 of eqn. (0) therefore is a power series in the coupling constant - but still it 
includes a non-perturbative effect in form of the non-analytical behavior in the temperature, 
which is of order g'^T log(T). 

One might argue, that these non-analytical terms are not significant because in order g^ 
other two-loop diagrams also give a contribution that we have neglected. However, as may be 
expected from our results, their non-analytical contribution involving log(T) is suppressed 
with respect to their leading order contribution by a factor 7, which may be translated into 
a factor g'^T. This implies, that our results to order g^ yield the dominant non-analyticity 
and are important for small temperatures (see figures). 

Moreover, this chain of arguments also supports the conclusion that terms of order 
g'^Tp^ /M'^\og{l/ g) cannot appear: Terms logarithmic in the coupling constant wear at least 
a coefficient g^ or g'^T'^/M'^, even if vertex corrections are introduced in the calculation of 
the self energy function. 
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Figs. H and |^ depict the self consistent solution of eq. (Q) for scalar bosons coupled to 
the fermions. For small temperatures, the deviation from the linear temperature dependence 
is very striking: Due to the Tlog(T) terms in the self-energy functions, 75(T) rises sharply 
at very small temperature (see comment in the next section). 

Before discussing this in detail, we turn to the results for the vector and pseudoscalar 
boson exchange. The self-consistent 7v due to massless vector boson exchange in Feynman 
gauge is plotted in figs. ^ and |^. The curves basically employ the same features as for the 
scalar bosons: A pronounced non-analyticity of the function jviT) at small temperatures. 

For the pseudoscalar boson our approach leads to a result which is surprising at the 
first glance: The self- consistent solution for jp is identically zero at all temperatures. This 
follows from the fact that the leading terms in the temperature dependence of eqs. (^) and 
(^) differ only by a factor of M/uo. Inspection of eq. (POj ) then reveals that these leading 
terms cancel, consequently /o = in the above equation ( ^41) and therefore 7p = 0. 

However, at a second glance this result is not surprising. The pseudoscalar coupling is 
a p-wave coupling, hence vanishes for zero momentum transfer. This is reflected in the fact 
that the coefficients /i and /2 start with power instead of possessing a constant piece. 
We therefore expect a self-consistent 7p which also starts with this power - clearly in 
contradiction to our ansatz spectral function. 

Naturally it is possible to repeat the calculation with a corresponding ansatz for 7 - 
but this would overstress the goals of the present paper. We therefore use this pseudoscalar 
result merely as a consistency argument in favor of our approximations. 

Also another feature of our results for 75 and 7y may be seen in the figures. Clearly, the 
momentum dependence of our self-consistent solution is small up to quite high temperatures. 
In view of eq. ( pS]) this means |t'o|>,|'^i| for scalar and vector boson exchange. We therefore 
conclude, that the ansatz of a momentum independent 75 and 7^/ is very well justified. 



A. Analytical approximations for the damping rate 



It is not necessary to repeat the analytical expressions for the self energy pieces in order 
to give the three functions /o - /2 in closed form. However, it is instructive to extract the 
dominant terms in the calculation of 7. 

To this end we perform an expansion of the solution of (|5^ ) in powers of the coupling 
constant. Only even powers occur, hence the expansion parameter is 



a 
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The contributions of order T and of order T log(T) then are, to second order in a 
2 p2 



75 ~ aT \ 1 



3M2 



+ 



a — log 

TT 



11 p2 



+ 



1061 



12 M2 1440 



12 



2^ 
3M2 



p 



15 



p6 rp2 



(55) 



(56) 
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where p = \p\. Naturally this last expansion is only good in the weak coupling limit. 
Performing the same expansion for the vector boson exchange yields 



-fv^aT 1 - 



2 p2 s pi 
SAP ^ 15 M4 

71 \ M\ 3 h'P 180 M' 



2 p2 8 
" 3M2 ^ 15M4 



p6 rj,2 



(57) 



It is a crucial aspect of these results, that they are analytical functions of the coupling 
constant, i.e., in contrast to other authors we do not find contributions of order \og{g'^) 



During the discussion of the self-consistency criterion ( pQ]) we have already stated, that 
the vector boson result is gauge independent in the limit p — 0. A vertex correction is 
needed to ensure gauge invariance for the momentum dependent parts (which in principle 
also require a completely different calculation scheme). However, this will at most give a 
regular modification of order a^p^/M^, logarithmic terms may appear only in even higher 
order a. 

The contribution of order T log(T), apart from being non-analytical around T = 0, has 
the effect of destroying the linear relationship between 7 and T: Neither for high tempera- 
tures nor for low temperatures can one approximate the self-consistent 7 by a linear function. 
This functional piece constitutes the leading non-analyticity observed numerically in figures 

(I) - m- 



B. Where does the log go ? 

A striking feature of our result is, that for a slow massive fermion there are no contri- 
butions to the damping rate of order g'^T \og{l / g). This seems to contradict some existing 
papers on the subject of the fermion damping rate, eg. refs. whereas others clearly 

denounce the idea of a logarithmic correction [^,0. We therefore have to set our result in 
relation to such modifications. 

Since our study was to a large extent motivated by ref. |@], we take the liberty to quote 
their eq. (15), with e replaced by g: 

l^hi^ /irarctan - (58) 



277^ j k \ 7 





g\ fdk g\ (gT\ 

7 

~ T log ( - J according to ref. [8] . (60) 



47r \g 
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Before we discuss the difference of tliis approacli to ours, let us first point out that also in 
our scheme appears a similar expression: In the calculation of Tip as well as in Ti/ appear 



terms of the type 



„ T / 2kLj ± 2pk - -f^ 

ak — arctan — — 

k \ 2'y[u! ± k) 



-k/T 



(61) 



They are due to the combination of angular integral factors with the leading term in 
the expansion of the Bose-Einstein distribution function (|50D . Restricting the parameters 
contained here to the approximations p = 0, A; <^ cj, 7^ <^ 2kuj, one indeed obtains the same 



integrand as in eq 



Two questions arise now: First of all it is important to ask, whether the "standard" result 
from eq. (|60|) is correct in itself, i.e., whether the final expression really is an approximation 
to the first expression. To answer this in detail, we note that one may indeed evaluate the 
integral in eq. (|58| ) in closed form. It is obtained as the Lerch transcendent, a generalization 
of the polylogarithm function and not expressible as a finite series of elementary functions. 

However, its numerical evaluation is straightforward, and we find that with reasonable 
accuracy it is indeed approximated by the expression in (|59|) - but not by the final expression 



(pDI). Instead, it follows from that 




log 



'Or. 



log ( - 

9 



-V 



(62) 



where e = 2.71828 .... As can be seen from fig. 8, this makes a tremendous difference: Even 
for coupling parameters as small as in QED, the step from (|59D to ( |60D is not justified, and 
the difference between the two expressions is not of higher order. 

The second question, equally relevant for the present paper, is the consistency of the 
first line in eq. (^). To this end we notice that in our approach there is no restriction to 
the momentum k: The fc-integration is automatically cut off at a scale k T due to the 
exponential function, see the expansion in eq. (|50|). 

Even more, as we have found the zero temperature pieces (involving 0-functions) of the 
distribution functions are crucial in replacing the non-analyticity in the coupling constant 
by a non-analyticity in the temperature. The answer to the second question therefore is, 
that the first line in eq. (^8]) (and similar expressions) may be consistent in the framework of 
hard thermal loop calculations, but is not a good approximation to the fully self-consistent 
damping rate problem. 

We may summarize the difference of the two approaches by symbolically writing down 
two different approximation schemes for the distribution functions and the corresponding 
result for the imaginary part of the self energy function: 



hard thermal loop resummation: 

nB{k) + npioj + k) ^ ^e{gT-k) 
UBi-k) + np{uj ~ k) ^ -lQ{gT + k) 

this work: 

ueik) + upiu + k) ^ (f + i) exp(-A;/T) 
nB{-k) + upiuj -k) ^ - (f + exp(-A;/T) 

-e(cu - k) 



— Tlog 

An \ 7 
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Naturally, in our approximation scheme one may no longer use -C but has to perform the 
full calculation as presented in the foregoing sections. We therefore conclude at this point, 
that the absence of terms of order g'^T\og{l/g) in our analytical approximations cannot be 
used as a reliable argument against our work. 



VII. CONCLUSION 

With the present paper we have reached several goals. First of all, we have presented 
a straightforward application of the method of generalized free fields. This automatically 
accounts for the undeniable mathematical fact that at nonzero temperature a perturbative 
expansion in terms of quasi-particle states with infinite lifetime is ill-defined [13-|l5[]- 



As we have shown, this method is also practical in the sense that it removes any unphys- 
ical infrared divergence from a perturbative expansion. We may put this into the following 
form: In reality a fermion in a heat bath is always subject to some Brownian motion, hence 
will never be in the same state long enough to emit a very soft boson. 

The central ingredient of the method are (more or less) continuous spectral functions. 
These functions are not a priori fixed in a theory involving generalized free fields ||15[, one 
may determine them from experiment or use physically reasonable self-consistency schemes. 
One such scheme we have presented here: In assigning a constant spectral width to fermions, 
we go one step beyond the introduction of temperature dependent ( "thermal" ) masses. 

This approximation is reasonable for systems which are dominated by the low-momentum 
sector of the interactions, i.e., those which involve massless gauge bosons. Naturally, to 
determine transport properties of plasmas in true non-equilibrium states, all the details of 
the spectral function have to be known and the approximation of a constant damping rate 



will break down 121 



For the massless bosons we used a "free" spectral function, without even a constant 
spectral width parameter. Although, in the spirit of our previous remarks, one may expect 
this to lead to inconsistencies in higher orders of perturbation theory we have made this 
choice for two particular reasons. First of all our calculation clearly shows, that and exactly 
how this worst-case infrared divergent interaction (which includes the Coulomb singularity) 
is regularized by the use of a nontrivial fermion spectral function. Secondly, any more 
realistic spectral functions for the bosons will only weaken the infrared behaviour of this 
spectral function, therefore leads to a smaller self-consistent damping rate. 

Numerical calculations show, that using such damped bosonic interactions in the self- 
consistent calculation is equivalent to inserting the sum of the fermionic 7 and the bosonic 7^ 
on the right side of eq. (0), i.e., in the calculated imaginary parts of the self energy function. 
Such a substitution does not affect the lowest order thermal contribution as obtained in sect. 
V. of this work. Thus it must be concluded, that the use of an effective boson propagator 
does not change the leading order result for the fermionic spectral width. 

For our approximate fermion spectral function we have achieved the summation of all 
nested Fock diagrams to arbitrary order. Apart from vertex corrections, such a summation 
constitutes the solution of the full one-body problem. As we have pointed out, such vertex 
corrections would lead to momentum dependent corrections to the damping rate - and 
therefore would not affect the momentum independent parts of our results. 
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In the weak coupling limit these momentum independent parts are analytical functions of 
the coupling constant, which sets them apart from some published results obtained with the 
method of hard thermal loops. We therefore made a careful investigation of the reason for 
logarithmic corrections within this method, and found them to be due to a certain approxi- 
mation made to the Bose-Einstein and Fermi-Dirac distribution function. A comparison to 
the approximation used in the present work is summarized in eq. (|63D . 

We find that logarithmic (i.e., \og{l/g)) corrections may appear only in higher orders of 
the self-consistent calculation, i.e., when including the self energy function to order 7^ log(7). 
As we have demonstrated by our use of a "free" boson propagator, they are certainly not 
due to the absence of screening on boson interaction lines. 

Another crucial difference of the calculations presented here to those established in the 
literature is the use of entirely causal propagators. Our spectral functions obey the locality 
axiom of quantum field theory (sect. II. B and D), whereas numerical studies show that this 
is not the case for those generally used in "hot gauge theory" One may argue, that 



this violation is of minor importance for the problem of the damping rate - but in our view 
the interaction of a fermion with its own radiation field should be strictly causal. 

Instead of a non-analyticity in the coupling constant, our results exhibit a non-analyticity 
in the temperature parameter arising from contributions of order Tlog(T). This finding is 
in accordance with the symmetry of space and time: In the presence of temperature the 
Poincare symmetry is broken to SO(3)xT4, and the symmetry restoration with temperature 
T — > is expected to be singular. 

We must stress the fact, that the calculation of the real part of the self energy function 
by a dispersion integral was suppressed in the present paper. According to the arguments 
we presented, the dispersion integral may still lead to non-analytical contributions in the 
coupling constant. However, these we expect to be of order exp{—l/g^) due to the removal 
of the Landau ghost poles 0] , and therefore they are varying slowly around the zero point 



of the coupling constant. 

A final conclusion to be drawn from our paper is, that it is certainly dangerous to 
present non-analytical results, like e.g. of order 5f^Tlog(T), without elaborate cross-checks. 
We have solved this problem by performing numerical calculations along with the asymptotic 
expansion of loop integrals, hence each figure presented here contains results obtained by 
two independent methods. Furthermore, every expansion made was carried out to one order 
higher than was finally used. 
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APPENDIX A: ANGULAR INTEGRATION 

In the following we give the explicit form of the angular integrals over the approximate 
fermion spectral function, which occur in the calculation of (^1]) and (^2|). rj in this appendix 
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resembles the cosine of the angle between the two momenta k (internal) and p (external), 
and we will henceforth write k = \k\, p — \p\. 



h 



1 

{k±po) = J 



dri ■ 



-1 



ci — 2pkf]y + cl 



i.pk{k ± po)7 



arctan 



{k±pof-t^-^' 
2{k±poh 



<^p-k 



and 



l2{k±po) = Jdv 



-1 
2kp 



Vcs + 2pkri 
Cl — 2pkr]y + c| 



(Al) 



(A2) 



Ii{k±po) + 



Sp'^k'^ 



log 



'{t^ + j^-{k±porf + A'y\k±pof 



r 



with 



C2 = 2(A; ±po)7 
C3 = a;^ + A;^ 



and boundaries of the integration defined as 



to. 



'l^, = u'±2pk + k' . 



(A3) 



(A4) 



An important aspect of the above analytical results is to select the proper Riemann sheet 
for the arctan-functions: /i and I2 are continuous functions in each variable. 
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FIGURES 



FIG. 1. Imaginary part F^^^ of the self energy function. 
Temperature T = 0, momentum p = 0, factor g'^/{47r^) set to one. 
Different constant values of 7: Thin line: 7 = 0, eqn. (^5]); thick line: 7 = 0.2M 

FIG. 2. Imaginary part F^ac: Numerical calculation vs. analytical approximation. 
Top panel: vector piece F^^^; bottom panel: scalar piece F^£ 

Temperature T = 0, energy po = u, momentum p = 0.5M, factor g'^/Air^ set to one. 

Continuous thin line: Full numerical calculation. 

Dash-Dotted thick line: Expansion to order 7. 
Dashed thick line: Expansion to order 7^. 
Continuous thick line: Expansion to order 7^. 

Dotted thin line: Expansion to order 7^, but without the term of order 7^. 

FIG. 3. Imaginary part F^: Numerical calculation vs. analytical approximation. 
Top panel: vector piece F^ ; bottom panel: scalar piece F^^'. 

Energy po = u, momentum p = 0.5M, factor g^/{47T^) set to one. 
Temperatures are T = 0.05M (bottom lines) and T = O.IM (top lines). 
Thin continuous lines: Full numerical calculation. 
Dash-dotted thick lines: Asymptotic expansion to order 7 
Dashed thick lines: Asymptotic expansion to order 7^ 

Thick continuous lines: Expansion to order 7^. 

FIG. 4. Self-consistent 75 /M as function of T/M 
Strong coupling a = 1, compared to 7 = T (thin line). 
Two different momenta: p = (thick continuous) p = 0.5M (thick dashed) 



FIG. 5. (75 - aT)/M as function of T/M 
Top panel: strong coupling a = 1; bottom panel: weak coupling a = 0.1 

Continuous lines: p = 0; dashed lines: p = 0.5M 

FIG. 6. Self-consistent 7y/M as function of T/M 
Strong coupling a = 1, compared to 7 = T (thin line). 
Two different momenta: p = (thick continuous) p = 0.5M (thick dashed) 

FIG. 7. (7y - aT)/M as function of T/M 
Top panel: strong coupling a = 1; bottom panel: weak coupling a = 0.1 

Continuous lines: p = 0; dashed hues: p = 0.5M 
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FIG. 8. Comparison of analytical approximations to eq. (58). 
Continuous thick line: numerical calculation of (ISSf) . 
Upper dashed line: approximation according to ([591). 
Lower dashed line: approximation according to (|60D. 



Continuous thin line: approximation according to (|62D. 
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FIG. 1. Imaginary part T^^^ of the self energy function. 
Temperature T = 0, momentum p = 0, factor g'^ / (Att^) set to one. 
DilFerent constant values of 7: Thin line: 7 = 0, eqn. (45); thick line: 7 = 0.2M 
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FIG. 2. Imaginary part F^ac: Numerical calculation vs. analytical approximation. 
Top panel: vector piece ; bottom panel: scalar piece F^^ 

11 1 vac ' 1 1 vac 

Temperature T = 0, energy po = cu, momentum p = 0.5Af, factor g'^/Air^ set to one. 

Continuous thin line: Full numerical calculation. 

Dash-Dotted thick line: Expansion to order 7. 
Dashed thick line: Expansion to order 7^. 
Continuous thick line: Expansion to order 7"^. 

Dotted thin line: Expansion to order 7"^, but without the term of order 7^. 
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0.05 0.1 0.15 0.2 0.25 

FIG. 3. Imaginary part F^: Numerical calculation vs. analytical approximation. 
Top panel: vector piece F^'; bottom panel: scalar piece F^^'. 

Energy po = cu, momentum p = 0.5Af, factor (/^/(Ftt"^) set to one. 
Temperatures are T = 0.05Af (bottom lines) and T = 0.1 Af (top lines). 
Thin continuous lines: Full numerical calculation. 
Dash-dotted thick lines: Asymptotic expansion to order 7 
Dashed thick lines: Asymptotic expansion to order 7^ 
Thick continuous lines: Expansion to order 7"^. 
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FIG. 4. Self-consistent ^js/M as function of T/M 
Strong coupling a = 1, compared to 7 = T (thin line). 

Two different momenta: p = (thick continuous) p = 0.5M (thick dashed) 
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FIG. 5. (75 - aT)lM as function of T/M 
Top panel: strong coupling a = 1; bottom panel: weak coupling a = 0.1 

Continuous lines: p = 0; dashed lines: p = 0.5Af 
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FIG. 6. Self-consistent -jv/M as function of T/M 
Strong coupling a = 1, compared to 7 = T (thin line). 

Two different momenta: p = (thick continuous) p = 0.5M (thick dashed) 
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FIG. 7. (7y - aT)/M as function of T/M 
Top panel: strong coupling a = 1; bottom panel: weak coupling a = 0.1 

Continuous lines: p = 0; dashed lines: p = 0.5M 
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FIG. 8. Comparison of analytical approximations to eq. (58). 
Continuous thick line: numerical calculation of (58). 
Upper dashed line: approximation according to (59). 
Lower dashed line: approximation according to (60). 
Continuous thin line: approximation according to (62). 



9 



